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1. Introduction 

Very recently, the Belle collaboration observed two narrow peaks, which are named 
Zb(10610) and ^{,(10650), in the 7r±T(n5) {n = 1,2,3) and 7c^h{mP) (m = 1,2) 
invariant mass spectra in the hidden-bottom decay channels of T{5S) The measured 
masses and widths of the two structures are, 

Mz,(io6io) = 10608.4 ± 2.0 MeV, T = 15.6 ± 2.5 MeV 
Mz,(io650) = 10653.2 ± 1.5 MeV, T = UA ± 3.2 MeV. 

Analysis favors quantum numbers of I^{J^) = l^(l^) for both states. The ^(,(10610) 
and Zb( 10650) are both charged bottomonium-like resonances and the masses are very 
close to the thresholds of the open bottom channels 5*5(10604.6 MeV) and B*B* 
(10650.2 MeV), so the molecular states of S-wave B*B and B*B* assignment are 
suggested by Belle collaboration. 

In the hadron level. Nils A. Tornqvist investigated the deuteron-like meson-meson 
bound states by meson exchange model [2]. The study shows that the energy of isoscalars 
BB* with J^^ = 0-+,l++, B*B* with J^^ = 0++, 0"+, 1+", 2++ are about 50 MeV 
below the corresponding BB* and B*B* thresholds. No bound state, however, appears 
for isovectors. Recently, by taking the pseudoscalar, scalar and vector mesons exchange 
into account in the framework of the meson exchange model, Liu et al. found that the 
loosely bound states probably exists in S-wave BB* [3l H]. Very recently. Sun et al. 
believe that the Zf,(10610), ^^(10650) are respectively B*B and B*B* molecular state 
after considering S-wave and D-wave mixing [5, 6j. 

In the quark level, by solving the resonating group method equation, Liu et al. [7] 
also investigated the system composed of {q = u,d,s) 

by two chiral quark models in which the pseudoscalar, scalar and vector mesons 
exchange are taken. The isoscalars BB, BB*{C = +), B*B*{J = 2) favor molecular 
states. Bondar et al. also discussed the heavy quark spin structure of the Zb(10610) 
and Zf,(10650) assuming that these are molecular state B*B and B*B* [8]. By 
considering the contribution from the intermediate Zf,(10610) and Zb(10650) states to 
the T(5S') — i- T(2S')7r+7r^ decay process, the anomalous T(2S')7r+7r^ production near 
the peak of T(5S') at i/i = 10.87 GeV [9j, observed by Belle collaboration, can be 
explained naturally [10]. The possibility of Zb(10610) and Zb(10650) being tetraquark 
states are discussed by the authors of Ref [HI |12] . The authors of Ref [131 E] • think the 
tetraquark and molecular structure both can interpret the Zb(10610) and ^(,(10650) in 
the QCD sum rule calculation. Further theoretical efforts concern the decay and mass 
of the Zb(10610) and Zfe(10650) states discussed in Refs. [l5l[T6l[n]. 

Inspired by the new states Zb(10610) and Zf,(10650) reported by Belle 
collaboration P and the related work, a systematical study of the possible S-wave 
BB, BB* and B*B* states is performed in this work. Here the B and B stand for 
{B^,B^,B^) and {B^,B^,B^^) triplets, respectively. It is worthwhile to investigate the 
intrinsic structure of the ^^(10610), ^^(10650) and other possible exotic states with h,h 
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quarks, especially in view of the great potential of finding new particles at Belle, BaBar, 
LHC and other collaborations. 

To study the mass spectrum of above possible exotic states, two types of chiral quark 
models (ChQM) [I8] are employed in this work. The numerical method, which is able 
to provide almost exact solutions, is very important in the study of few-body systems. 
Here, a high precision numerical method for few body system, which is different from 
the methods used in the previous work by other researchers, the Gaussian Expansion 
Method (GEM) is used. The detail of GEM can be found in Refs. ^ ED]. 

The paper is organized as follows. In the next section we introduce the Hamiltonian 
of the chiral quark models. Section |3] is devoted to discuss the wave function of possible 
molecular states BB, BB* and B*B*. In Section HI we present and analyze the results 
obtained in our calculation. Finally, the summary of the present work is given in the 
last section. 



2. The chiral constituent quark model 

In the ChQM, the Hamiltonian usually includes Goldstone-boson exchange in addition to 
color confinement and one-gluon-exchange (OGE). The chiral partner, cr- meson, is also 
usually introduced, although its existence is still in controversy [21j. The Hamiltonian 
of the ChQM used here is given as follows, 

H=j:(m+^]-TcM+ i: {vs + + + (1) 

i=i \ -^rni/ j>i=i 
where x = 7r,i^, 77, Tcm is the kinetic energy operator of the center-of-mass motion of 
whole system. 

The linear confining potential, which is suggested by lattice QCD calculation of qq 
system, can be written as 

Vg = X'r-X'j{-a,r,,-A). (2) 

For one-gluon-exchange, the potential takes the form 

1 277 



o ■ ^j) ^(^ij) 



(3) 



where, cr, A are the SU(2) Pauli matrices and the SU(3) Gell-Mann matrices, 
respectively. The A should be replaced by —A* for the antiquark. In the non-relativistic 
quark model, the delta function 6{rij) should be regularized |22], because of the finite 
size of the constituent quark. The regulation is flavor dependent and reads [HI |23] 

^i^^j) = 7 ^ e-^-/^^°('^), (4) 

where ro(/i) = ro//i and /i is the reduced mass of quark-quark (or antiquark) system. 
The wide energy covered from light to heavy quark requires an effective scale-dependent 
strong coupling constant in Eq. ([3]) that cannot be obtained from the usual one-loop 
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expression of the running coupling constant because it diverges when Q — )• Aqcd- Hence 
an effective scale-dependent strong coupling constant [T8] is taken as 

"0 



In 



^^0 



(5) 



where no and Aq are the free parameters. 

For the mesons exchange, potential takes the form 

.2 

-Hi{m^,A^,rij){(Ti 



ml 



12m,mi 



a=l 



Vif = C{g,h,AK,mK) 



m 



K 



12m,m. 



(6) 
(7) 



a=4 



C{gch,A.^,mrj)- 



12mimj 
cos^p(Af ■ Af) -sin^p(A°- A' 



-C{gch,A^,m, 



Hi{m, A, r) 
H2{m,A,r) 
C{gch,A,m) 



Y(mr) 



A^ 

— F(Ar) 



Y(mr) 



9ch 



A' 



—Y(Ar) 
m 



-m 



(8) 
(9) 

(10) 

(11) 
(12) 



Att A^ — m? 

where the a exchange only occurs between the lightest quarks {u- or d-quark) which is 
different from Ref. [18] due to its non-strange nature. The strange scalar meson (with 
large mass) exchange is not taken into account in the present work because of its small 



effect. We adopt A° 



/ due to the normalization of SU(3) matrix. Y{x) is the 



standard Yukawa function defined by Y{x) = e~^/x and the rest symbols have their 
usual meaning. The chiral coupling constant g^h is determined from the ttNN coupling 
constant through 



9l^ 



3x2 „2 



9nNN ^u,d 

Att 



m 



(13) 



N 



and flavor SU (3) symmetry is assumed. The tensor term and the spin-orbital term have 
been omitted in the potentials since we consider only S-wave states. 

The above model is denoted as ChQMl. To testing the effect of cr-exchange between 
the lightest and strange quark or strange quark pairs, and long-range color screening 
on the binding energy of the molecular states, the Salamanca version of the chiral 
quark model [18], which is referred as ChQM2, is also employed here. The screened 
confinement interaction in this model is 

Vg = Xt ■ \] {-a,(l - e-'^--) + A}, (14) 

where /ic is a color screening parameter. The other potentials are the same as the above 
with the exception that the cr-meson is exchanged between all the light quarks, m, s. 
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3. Wave function 

The total wave function of multi-quark system can be written as, 

^i±= lOlv^' ^JJ., (15) 



where |^), \r])\ \x)s^ \^) Lt I'spresent color singlet, isospin with /, spin with S and 
spacial wave function with angular momentum Lt, respectively. 

All possible molecule structures composed of S-wave B and B, which stand for 
[B'^jB^jB^) and {B~,B^,B^) triplets, respectively, are investigated in this work. 
According to the total isospin, the PP (pseudoscalar meson) and VV (vector meson) 
flavor wave functions are listed in Table [H Another possible molecule structure for 
four-quark system is bottomonium-|-light meson. In this case, there is no interaction 
between colorless bottomonium and light meson (color dependent interaction is zero 
between two colorless cluster if no exchange term exists and there is also no Goldstone- 
boson exchange between heavy and light quarks), so no bound state can be formed in 
this case. Therefore, we do not take into account of this case in the present work. 

Obviously the components PV and VP do not have definite C parity, one can get C 
parity = ± by [PV ±C (PV)]/ V2 and [VP±C {VP)]/V2 for the neutral states |2] such 
as (I, /^) = (1,0), (0,0) shown in Table [H (All the orbital angular momenta are set to 
zero because we concentrate on ground states). Hence, the coefficient C = ±1 represent 
C-even and -odd parity respectively, which are different from that of Ref. [U IH |2^ . 
since we use normal convention of PDG p5] i.e. B^ = db and B^ = db. One can 
easy find that G parity = =F are corresponding to C parity = ± for these states with 
(I, J^) = (1,0). However, there is no interaction depending on C and G parity in 
the model Hamiltonian Eq.([T]), so these two states with ± C and G parity must be 
degenerate in our calculation. The two states separated by comma in each row of Table 
[Uwill be coupled in the calculation. 

The spatial structures of molecular states are pictured in Fig. [1] The relative 
coordinates are defined as following, 

r = ri - r2, R = rg - r4, (16) 

P = , , , (17) 

nil + ^2 ^3 + fn^ 

and the coordinate of the mass-center is 

4 4 

R-cm = ^i^il H "^i' (IS) 
i=l i=l 

where is the mass of the ith quark. 

Then the outer product of space and spin wave functions is 

Xs2nis2]j2M2]ji2Mu'^%{p)]jJ,- (19) 
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Table 1. The flavor wave functions of the BB, BB* systems. "C" is the charge parity. 



Isospin 


BB 


BB* 




B+B° 


B*+B°, B+Bf 




B°B° 


B*^B°, B°Bf 




B^B'^ 


BfB°, B^^B*^ 




B°B- 


BfB-, B^^B*- 




B+B" 


B*+B", B+B*° 



1=1 -^{B+B- - B°B°) ^[{B*+B- - B*°B°) + C{B*-B+ - B*°B°)], 

1[{B+B*- - B°B*°) + C {B-B*+ - B^B*°)] 
B"B- B*°B-, B^B*- 



1=0(0 


7l(^' 


B- +B°B°) 


i[(B*+B-H 


^ B*°B°) + C{B*-B+ ^ 


-B*^B°)], 






^[{B+B*- - 


^B"B*°) + C{B-B*+ H 


-B°B*°)] 


l=0(s) 








hC4*0i?0), j^iB^^Bf 


+ B°Bf) 




Figure 1. The relative coordinate for a meson and antimeson system, solid and hoUow 
circles represent quarks and antiquarks, respectively. 



In GEM, three relative motion wave functions are written as, 

0L(r) = E CnN^ir'e-'-^^Ui) (20) 

n=l 

^ max 

^lmCR) = E c^iV^Li?''e-'^-«ViM(R) (21) 

N=l 

^%{P) = E c^N^ppPe---''%,{p) (22) 



a=l 



Gaussian size parameters are taken as geometric progression 

1 

1 



^^r^ = r,a-\a=l^y- ' (23) 



The expression of (jy, RN,<^a, Pa in Eqs. (12T]) - ( 122|) are similar to Eq. (123|) . 

The physical state must be in color singlet, which can be constructed in two ways: 
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color- singlet and color octet, 

l^l) = |ll2® I34) 



16) = |8i2 ® 834) . 



(24) 



The state in color octet channel is called hidden color states by analogy to states 
which appear in the nucleon-nucleon problem [26] . 

The total spin of PP and PV and VP system is only 0, 1 respectively. However, 
the VV system can coupling to total spin 0, 1 and 2. 

4. Numerical results and discussion 

Solving the Schrodinger equation 



with Rayleigh-Ritz variational principle, the energies of normal mesons, BB, BB* and 
B*B* systems can be obtained by using different total wave functions, respectively. 

To determine if the BB, BB* and B*B* systems are bound or not, the threshold of 
the system should be fixed. Clearly the threshold is governed by two corresponding 
meson masses. So one believes that a good fit of meson spectra, with the same 
parameters used in four-quark calculations, must be the most important criterium|23[ 
[28l [271 l29l Eni [31]. Of course, there is another possible threshold for four-quark system, 
bottomonium+light meson, e.g., T(IS') + p for IJ^ = 11"^ channel. Generally the 
threshold in this case is lower, so the bound state B^*^B^*^ will become a resonance in 
this channel. Since the transition from B^*^B^*'^ to bottomonium+light meson involves 
string rearrangement, we leave this for the future work. 

In GEM, the calculated results of normal meson spectra (or the spectra of BB, 
BB* and B*B* systems) are converged with the number of gaussians Umax — 7 {rimax — 
7, Njnax = 7, amax = 12), and the size parameter r„ (r„, Rjq, p^) running from 0.1 to 
2 (2, 2, 6) fm. The convergence properties of the energies have been discussed in detail 
in Ref.[20]. The parameters and the normal meson spectra in two types of ChQM are 
listed in Table [2] and [3l respectively. 

By solving the equation f l25l) . the energy of the BB, BB* and B*B* systems can be 
obtained. If the binding energy, /S.E = Mgystem — — (g = u, d, s), is negative, 
then the system would be bound. According to the Table. [HI the thresholds of possible 
molecular states BB,BB*,B*B* are easily listed in Table [H 

The color-singlet, color-octet channel, and channel coupling calculation of 
BB, BB*, B*B* systems are done in the two types of ChQM, and the results are presented 
in Table [MZl in which the results are denoted by "1®!", "8® 8" and "channel coupling", 
respectively. Due to the strong interaction is invariant under the rotation of isospin, the 
different states corresponding to the different components of isospin I are degenerate, 
so we present the results for each total isospin /. From Table [5}[Zl we can see two models 
give very similar results. 

In the S-wave PP system, the quantum numbers are O"*". Apart from scalar 
meson a, the pseudoscalar mesons e.g. tt, K, rj cannot be exchange in BB system because 



{H - E) vl>^/x = 



(25) 
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Table 2. Parameters of two quark models. The masses of tt, 77 in Eqs.(|6l)-([8|) are got 
from the experimental data which are = 0.7 fm"-'^, m,, = 2.77 fm^^, respectively; 
TOo-,A„, A^, Op are taken the same values as Ref. jl8j. namely — 3.42 fm~^, 
= = 4.2 fm-i, A^ = 5.2 fm-\ 0p = -15°, g2^/47r=0.54 







ChQMl 


ChQM2 [18] 


Quark masses 


mu,d (MeV) 


313 


313 




ms(MeV) 


525 


555 




mc(MeV) 


1731 


1752 




mb{MeV) 


5100 


5100 


Confinement 


ac(MeV fm-i) 


160 


430 




A (MeV) 


-131.1 


181.1 




/Xc(fm"^) 




0.7 


OGE 


ao 


2.65 


2.118 




ro(MeV fm) 


28.17 


28.17 




Mo (MeV) 


36.976 


36.976 




Ao(fm) 


0.075 


0.113 



Table 3. Numerical results of normal meson spectrum(unit: MeV) in the two quark 
models. The experimental data marked "Exp." takes from the latest Particle Data 
Group |25], and the ground state of bottom ?7f,(ls) is observed by BABAR Collaboration 
from the radiative transition T(3iS') — ^ 777b |32] 



Meson 


ChQMl 


ChQM2 


Exp. 


Meson 


ChQMl 


ChQM2 


Exp. 


TT 


140.1 


153.2 


139.57±0.00035 


Ds 


1966.6 


1991.8 


1968.49±0.34 


K 


496.2 


484.9 


493.677±0.016 


Dl 


2091.1 


2094.1 


2112. 3±0. 5 


p(770) 


775.3 


773.1 


775.49±0.34 


i?± 


5284.7 


5277.9 


5279.15±0.31 


K*{892) 


917.9 


907.7 


896.00±0.25 




5284.7 


5277.9 


5279.53±0.33 


w(782) 


703.7 


696.6 


782.65±0.12 


B* 


5324.3 


5318.8 


5325. 1±0. 5 


(/)(1020) 


1016.8 


1011.9 


1019.422±0.02 


B's 


5360.6 


5355.8 


5366. 3±0. 6 


'7c(ls) 


2995.7 


2999.8 


2980. 3±1. 2 


Bt 


5403.6 


5400.5 


5412.8±1.3 




3097.6 


3096.7 


3096.916±0.011 




9384.6 


9467.9 


9388.9;^;3(stat) 




1882.2 


1898.4 


1864.84±0.17 


T(ls) 


9462.4 


9504.7 


9460.30±0.26 


D* 


2000.2 


2017.3 


2006.97±0.19 











of the parity conservation. One can find in the Table |5] that the a meson exchange do 
not contribute enough attraction to bind the BB system in color-singlet channel in the 
ChQMl. Due to the a-exchange also occurs between ss pair and us or ds pairs, loosely 
bound states of BB with I = 0(s), or | are obtained in the ChQM2. Noteworthily, 
there is no one-gluon-exchange between two separate mesons just in this channel. More 
bound states are formed if we take the coupling of color-singlet and color-octet channels 
into account. Obviously, in the color octet channel, the color-magnetic terms of OGE 
between two separate colorful mesons contribute attraction to BB system, because of 
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Table 4. Threshold energy of BB,BB*,B*B* in two quark models. (unit: MeV). 



J 


I 


Ml Mo 

1 1 J. 1 vj. 2, 








0,1 




10569.4 


10555.8 










10721.2 


10711.6 




1 

2 


s 1 s 


10645.3 


10633.7 




0,1 


B°B*+ 


10609 


10596.7 


1 





B7B*+ 


10764.2 


10756.3 




1 

2 


B-B*°/B+B*- 


10684.9 


10674.6 






B*-B^/B*+B- 


10688.3 


10678.4 




0,1 




10648.6 


10637.6 


0,1,2 





B*-B*+ 


10807.2 


10801 




1 

2 


B*-B*°/B*+B*- 


10727.9 


10719.3 



the requirement of total color-singlet of the state. In this case, due to the masses of 
u, d quarks are much smaller than the mass of h quark, the cross matrix between color- 
singlet and -octet channels of the color-magnetic interaction, which is in proportional 
to l/(mjmj), is very large. So the energy of each BE system is depressed by it, which 
was discussed in detail in Ref . [231 ES] . 



Table 5. The binding energy (unit: MeV) of BB. The " 1 (g) 1" , "8 (g) 8" and "channel 
coupling" represent BB in color-singlet, color-octet and coupling of color-singlet and 
color-octet channel, respectively. 







ChQMl 






ChQM2 




Isospin 


l®l 


8® 8 


channel coupling 


1 ® 1 


8(g)8 


channel coupling 


1=^ 


0.5 


51.5 


-17.5 


-0.2 


97.7 


-2.6 


1=1 





0.1 


-72.6 





45.6 


-29.9 


I=0(/) 





0.1 


-72.6 





45.6 


-29.9 


I-O(s) 


0.3 


77.4 


0.2 


-0.7 


137.5 


-1.5 



For the BE* system, the calculation results are listed in Table [61 The a, vr, rj mesons 
can all be exchanged in such systems. Two bound states are both found in two quark 
models in color-singlet single channel. For 1 = 1 state, the binding energy are both 
about —1 MeV with regard to BB* threshold in two quark models, and the distance of 
each pair quarks in ChQMl are 

= 0.63 fm, = 0.63 fm, 

/(7£y = 2.1fm, /(^=2.1fm, 

= 2.02 fm, = 2.18 fm. 

The distance between quark and quark or antiquark belonging to different mesons are 
larger than that between quark and antiquark belonging to the same meson, the results 
show a clear molecule structure. It is reasonable to interpret the state ^^(10610) reported 
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by Belle collaboration as the molecular state BB* with I{J^'^) = l(l+~). For this state, 
there are also other thresholds, T{lS)p and hf,7r [15]. The energy of the state is a httle 
higher than these thresholds. Because of the different color structures the states BB* 
and T{lS)p, hf,7r have, the transition involves the color structure rearrangement, BB* 
may appear as a resonance in the T[lS)p and hbir channels |3l]. The calculation of the 
transition, which is out of the scope of the present work, is left for future work. In our 
calculation, the / = state without strange quark is also a bound state with binding 
energy —12.1 MeV in two quark models. In the ChQM2, in addition to the above two 
states, the I = and 0(s), also form a bound state for the a-exchange contributes to 
these channels. The color-singlet singlet and hidden color channel coupling leads to that 
all the states are bound. 



Table 6. The binding energy (unit: MeV) of BB* . 


ChQMl 

Isospin 1 8 (E)8 channel coupling 


ChQM2 

1 1 8 8 channel coupling 


I=i 0.4 -21.2 -89.8 
1+ 1=1 -1.3 -88.1 -164.3 
1=0(0 -12.1 -47.2 -122.6 
l=0(s) 0.35 14.8 -47.3 


-0.2 28.9 -40.9 
-1.1 -35.5 -107.4 
-12.1 2.0 -69.8 
-1.2 77.5 -3.9 


Table 7. The binding energy (unit: MeV) ofB*B*. 


ChQMl 

Isospin 1 (X) 1 8 8 channel coupling 


ChQM2 

1 (E) I 8 (8) 8 channel coupling 


I=i 0.4 -111.8 -172.9 
0+ 1=1 -3.4 -194.4 -266.2 
1=0(0 0.4 -104.9 -175.9 
l=0(s) 0.3 -58.9 -112.2 


-0.5 -56.6 -108.4 
-3.0 -132.3 -195.5 
0.4 -49.8 -115.3 
7.6 -37.7 


I=i 0.4 -71.4 -129.0 
1+ 1=1 -1.2 -135.7 -201.5 
1=0(0 0.4 -94 -160.1 
l=0(s) 0.3 -34.9 -86.0 


-0.3 -22.5 -73.9 
-0.9 -83.2 -144.6 
0.5 -44.9 -107.6 
-0.2 26.5 -21.2 


I=i 0.4 -2.6 -56.4 
2+ 1=1 0.3 -37.6 -97.0 
I=0(/) -11.1 -74.3 -133.0 
l=0(s) 0.3 8.4 -41.0 


36.1 -19.3 
0.4 -1.9 -62.6 
-11.0 -35.7 -95.1 
-1.1 61.1 -5.7 



The S-wave B*B* systems have quantum numbers J^'^ = O"*"^, l"*"", and 2+"*" for 
the neutral states. The tt, rj, a mesons can all be exchanged for 1=0(/) and 1=1, while 
only rj is exchangeable for l=0(s) and | states. The a interaction is always attractive 
between the lightest u, d quarks. According to Eq. (Q, the vr-exchange is attractive 
for the states with I{J^) = l(O^), l(l^), 0(2+) and makes these states are all bound 
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states, which are shown in Table [7] in the color-singlet single channel calculation. The 
B*B* with I{J^) = l(l^) has binding energy about -1 MeV in two quark models, and 
the distance of each pair quarks are 

= 0.64 fm, = 0.64 fm, 

^) = 2.11 fm, = 2.27 fm. 

The assignment of the newly observed state ^^(10650) to a molecular state B*B* with 
I{J^'^) = 1(1^ ) is favored. In the hidden color channel, almost all the states, except 
the one with IJ^ = 02+ and hidden strange, are become bound. Again the channel 
coupling between color-singlet and hidden color channels makes all the states bound. In 
the present work, the S-D mixing of is not taken into account, the mixing will be 
important for states with energy on the threshold. From the calculation of deuteron, 
we estimate the S-D mixing will increase the binding energy of Zf, about 2 MeV. 

5. Summary 

In the framework of chiral quark model, a systematical study of the mass spectra of 
BB, BB* and B*B* systems is performed. The states BB* and B*B* with quantum 
numbers I{J^^') = 1{V ) are shown to be bound, which are respectively good 
candidates for the charged bottomonium-like resonances ^^(10610) and ^^(10650) newly 
observed by Belle collaboration. The color-singlet single channel calculation also shows 
that the states BB* with J(J^^) = 0(1++), B*B* with J(J^^) = 1(0++), 0(2++) are 
bound states with a few MeV binding energy. 

Recently Belle collaboration reported their high precision measurement of 
bottomonium mass: M[T(5S')] =10.87 GeV [9]. If the molecular states B*B* with 
I{J^) = 1(0+) really exist, it could be observed in final state T(15')p at the Belle, BaBar, 
LHC and other collaborations. Due to the phase space limitation, the isoscalar states 
BB*{J^ = 1+) and B*B*{J^ = 2+) may be observed in decays of excited bottomonium 
which above the T(55'). 

The (T-exchange plays important role for binding the BB, BB* and B*B* with I = \ 
and 0(s) in the ChQM2. To search these molecular states in the future experiment will 
test the contribution of the cr-exchange in the chiral constituent quark model. 

The hidden color channel effect is complicated in multi-quark systems. Here we 
extend directly the quark model for colorless cluster to the colorful cluster in the study 
of BB, BB* and B*B* systems. In our calculation, the color-octet channel plays a 
dominate role in producing deeply bound states. If the quark-antiquark interaction in 
color singlet can be extended directly to color-octet by Casimir scaling [35], then the 
OGE interaction will be attractive between two color-octet cluster in some multi-quark 
systems. So it is inevitably to produce deeply bound states for hidden-bottom states 
because of the too small kinetic energy. More experimental data on bottomonium-like 
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resonances are needed to check the Casimir scahng, and cast doubt on the vahdity of 
applying the chiral constituent quark model to the hidden color channel directly. 
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